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Abstract 

The Casimir-Lifshitz force is calculated between two inhomogeneous composite slabs, each made 
of a homogeneous matrix with spherical metallic inclusions. The effective dielectric function of the 
slabs is calculated using several effective medium approximations and we compare the resulting 
forces as a function of slab separation and filling fraction. We show that the choice of effective 
medium approximation is critical in making precise comparisons between theory and experiment. 
The role that the spectral representation of the effective medium plays in making a Wick rotation 
to the complex frequency axis is also discussed. 
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I. INTRODUCTION 



The prediction by H. B. G. Casimir [1] that two neutral parallel plates made of a perfect 
conductor will attract each other due to zero point fluctuations has gone from a theoretical 
curiosity to a force that has been measured with very high precision and that has implica- 
tions in several fields. Knowledge of Casimir forces is needed to search for deviations from 
Newtonian gravitation at short distances ^2 4 ] and to detect chameleon fields |5j. In applied 
physics and engineering the Casimir force plays a role in micro and nano mechanical devices 



12|. 



The generalization of the Casimir force to arbitrary materials was developed by Lifshitz 

HQ 

. In the Lifshitz theory the force originates from fluctuating electromagnetic fields in 
the macroscopic bodies. These fluctuating fields are related to the dielectric function of the 
materials via the fluctuation-dissipation theorem. Although a closed form expression for the 
dispersive forces can be given only for parallel plates separated by a gap, Lifshitz theory is 
general since it includes both the retarded and non-retarded limits. For perfect conductors 
Casimir's original result is recovered. 

The dependence of the Lifshitz formula on the dielectric function and magnetic suscep- 
tibility [15) can be used to modify the behavior of the force. For example, its sign can 
be changed by a suitable choice of the dielectric functions of the slabs and of the medium 
between them 16]. Also, its magnitude can be chan ged using thermally activated structural 



)hase transformations of materials like AglnSbTe [17| or TiOi , external magnetic fields 
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20] or by changing the charge carrier concentration of the slabs 21]. 



Composite materials provide additional issues in the evaluation of the Casimir force. 
Typical composites are obtained by combining two or more different materials such that at 
certain wavelengths the composite behaves as an homogeneous material with an effective 



dielectric function. Composites such as V0 2 
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aerogels 



24|, composites with nickel 



particles 25j and most recently meta materials 26H30] have been shown to change the 
Casimir force. 

In general, a composite material is defined as a mixture of N different materials each with 
a physical property (such as a dielectric function) pj, each occupying a volume fraction fj 
such that YljLi fj = ^- The problem then is to relate the average or effective value p to the 
component values. In addition to dielectric function, other physical properties pj of interest 
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are the magnetic susceptibility, their Lame constants, the thermal conductivity etc. In this 
work we will be concerned only with the effective dielectric function e which in general is 
frequency and wave vector dependent. 

Understanding the effective properties of composites remains an important problem. Ex- 



amples where this is an important issue include magnetodielectric photonic crystals 341 ]. 



metamaterials 



351 1. surfaces and meta surfaces |36|. Composite materials are also important 
in coUoid chenhstry, and there ha, been recent interest in the optica! Q and St— 
properties of such materials as a result of advances in the preparation of colloidal crys- 



tals and c 
materials 



38| , and in the preparation of DNA-linked gold nanoparticle superlattice 



usters 

adj. 

How many ways are there to find the effective dielectric function of a composite system? 
In a recent paper by Prasad [3l| at least 12 different expressions or mixing formulas for 
calculating the effective properties are presented and compared with experimental results. 
There is no general answer as to which is the best effective medium approximation (EMA). 
As stated by Grandquist and Hunderi, the merits of the various theories can be judged only 



when comparing with experimental data 32]. Comparisons with more fundamental theories 



in which the properties of the constituent particles are explicitly included is also possible, 



as was applied for gold nanoparticles by Lazarides and Schatz [33 ]. 

In this paper we calculate the Casimir-Lifshitz force between two composite slabs made 
of a dielectric homogeneous matrix with spherical inclusions of Au. The effective dielectric 
function of the slabs is calculated using several of the effective medium approximations. 
Depending on the choice of mixing formula, the resulting Casimir force changes significantly. 

In the first section we present a brief summary of the most common effective medium 
approximations. Also we discuss the spectral representation and its importance in the proof 
of Kramers-Kronig relations needed to evaluate the effective dielectric along the imaginary 
frequency axis. Finally we present calculations of the Casimir force for different mixing 
formulas, showing the variations in force with the choice of mixing formula. 

II. EFFECTIVE MEDIUM APPROXIMATIONS 

In this section we review some of the existing effective medium approximations without 
going through their derivation. The system we consider is an homogeneous matrix made of 
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a material with a dielectric function (the host) which has inclusions of a second material 
with a dielectric function q (inclusions). The inclusions occupy a volume fraction / of the 
matrix, also known as the filling fraction. The dependence on frequency of the dielectric 
functions will be suppressed. Without loss of generality we will assume spherical inclusions 
of radius a. 



The simplest approach used by Wiener 40j was to take simple arithmetic and harmonic 
averages, that by definition will provide an upper and lower bound for e. This formula is 



< e < fe, t 



f)e h . 



(1) 



fe h + (1 - f)e 

More rigorous bounds for the effective dielectric function were found by Hashin and 



Strickman 



41] using variational procedures. Although the original derivation of these bounds 



was done for the magnetic susceptibility they are applicable to the dielectric function, electric 
conductivity and heat conductivities. The Hashin- Strickman bounds (H-S) are 



_j e ±l_ < 2 < f _i I 1 jj^± (2) 

h e h j{e i - e h ) + (1 - f)/3 ~ ~ ^ €i /(e h - e,) + //3" 1 ) 

The fundamental property of the H-S bounds is that any EMA has to yield a dielectric 
function that falls between these two bounds. It is important to notice that Eq. (1) makes 
no assumption on the shape of the inclusions. However the shape is important and the 
expression given in Eq. (2) is for spherical inclusions. 

The Maxwell- Gar nett approximation is the most commonly used EMA cited in the text 
books due to its simplicity and physically intuitive derivation 
function is 

Vol 



where a 



I- fa 

7- is the polarizability of the inclusion (divided by a 3 



40| . The predicted dielectric 

(3) 



The Maxwell-Garnett approximation (MGA) assumes that the separation between the 
inclusions is large so that they act as independent scatterers and that the polarization effect 
on one of the particles due to its neighbors is instantaneous. The first assumption limits the 
applicability of MGA to dilute systems with a low filling fraction /. This assumption breaks 



down as the filling 



been proposed 



32 



raction increases and corrections to the Maxwell-Garnett theory have 
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44j . The limitations and range of validity of the Maxwell-Garnett 



theory have been discussed extensively by Ruppin [45 1. 
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The Maxwell-Garnett theory does not predict percolation transitions. For example if 
the host material is an insulator and the spherical inclusions are metallic, as we approach 
the random packing filling fraction an insulator-conductor transition can be expected. An 



improvement to Maxwell-Garnett theory comes from Bruggeman |46| who treated the two 
components of the composite in a symmetrical fashion. In the Bruggeman approach the 
effective dielectric function is obtained by solving the equation: 

Both Eq.(|3]) and Eq. (j4j) are for spherical inclusions. However, they can be derived 
for other geometries. In the Appendix we present the Maxwell-Garnett and Bruggeman 
expressions for ellipsoids. Also, both expressions can be generalized to nonlocal dielectric 



functions wherein the dielectric function depends on frequency and wave vector. 47] 



As mentioned before, there are many EMA expressions and not all of them will be 
considered in this paper. Besides the Maxwell-Garnett and the Bruggeman EMA we selected 
Looyenga's approach 4^, which is in good agreement with experimental measurements of 
the effective optical properties of colloidal quantum dot systems [49]. Looyenga's model for 
spherical inclusions assumes an equation for the effective dielectric function of the form 

e=(fe} /3 + (l-f)el /3 )\ (5) 

again assuming spherical inclusions. The models so far presented are a small sample of many 
of the mixing formulas that exist. We emphasize that there is no generally accepted best 
effective medium theory so it is important to study a broad range of theories to see what 
sensitivity there is to the computer Casimir force. 



A. Spectral representation and Rotation to the Imaginary axis 



The effective dielectric function of a composite depends on the dielectric functions of 
the materials and on the geometry of the system, for example the shape of the inclusions. 
These two effects can be separated using the so-called spectral representation developed 
independently by Fuchs 50], Bergman and Milton 5l|, |52|]. By defining a function 



(6) 



that only has information on the material properties and a function G(L), the spectral 
function, that contains information on the geometry of the system, the effective dielectric 
function can be written 

where L e [0, 1] is the depolarization factor. The spectral function G(L) can be interpreted 



as the density of geometrical or shape resonances [53j]. For example, for the Maxwell- Garnett 
model for spherical inclusions one has that G(L) — 6(L — (1 — /)/3). 

Using the spectral representation one can derive the mixing formulas described in the 
previous section, including Winer bounds, Looyenga, Maxwell-Garnett and Bruggeman. The 
spectral function (dielectric constant) of the components can be determined experimentally 
from reflectivity measurements [5J], for example. 

There is a more fundamental consequence of the spectral representation. It can be proven 
that any effective dielectric function obtained using the spectral representation satisfies the 
Kramer-Kronig relations 55J • The importance of the Kramers- Kronig relations in Casimir 
force calculations is that the Lifshitz formalism relies on a rotation of the frequency to the 
complex plane, that is u — > i(. As shown in Ref. ([5y]) rotation of the dielectric function is 
obtained from the Kramers-Kronig relation as 



2 f , ue (oj) 

e(iO = 1 + - J d.-^. (8) 







B. Comparisons between the different models 

Consider a composite material that consists of a homogeneous host of S1O2 and spherical 
inclusions made of gold. The dielectric function of S'i0 2 is described by an oscillator model 

r/-\ 1 , Cuv G IR , s 



where Cir,uv are the a 
absorption frequencies 



Dsorption strengths in the infrared and ultraviolet and ojiruv the 



5t| . For the inclusions we consider spherical Au particles whose 
dielectric properties are obtained from tabulated data [58] with a low frequency extrapolation 
using a Drude model j^. In figure (1) we show the dielectric function calculated as a 
function of the filling fraction for the Wiener bounds (upper and lower bound), Maxell- 
Garnett, Bruggeman, and Looyenga expressions. The frequencies used ( = 0.02u;p (top) 



and C = 0.5w p (bottom) were chosen arbitrarily. For the top graph the effective dielectric 
function is plotted with a logarithmic scale. The large values of e are due to the pole 
at zero in the Drude dielectric function that is responsible for the DC conductivity. At 
this frequency the dielectric function of Au dominates. The Bruggeman model predicts a 
percolation transition at / = 0.33 that is evident from the inflection in the plot. For higher 
frequencies the dielectric function of Au decreases (and it should go asymptotically to 1 
with increasing frequency) and the dielectric function of the host medium is comparable to 
that of Au. In this case, the Maxwell-Garnett, Looyenga and Bruggeman models give very 
similar results as seen in the bottom graph. 



III. EFFECTIVE MEDIUM THEORIES AND CASIMIR FORCES 

To calculate the Casimir force we use the Lifshitz formula between two parallel slabs with 
local dielectric functions ei(uu) and €2(00) respectively, separated by a gap length L. The gap 
is filled with a dielectric function e 3 (w). Lifshitz result for the force per unit area is 

F -^Y. */ rf Q^ vJ-"r- ' (10) 

Z7r v=s>p J0 J0 e r 13 r 23 

where is the reflectivity between medium i and j for either p or s polarization, Q is the 
wave vector component along the plates, q = (/c and k 3 = ^/e 3 q 2 + Q 2 . The above expres- 
sion is evaluated along the imaginary frequency axis i(, and now the dielectric functions 
e i{K) (i — I5 2, 3) have to be considered. This expression is valid for hw > KbT, (Kb is the 
Boltzman constant). 

In our case we assume e 3 = 1 and the dielectric functions for both slabs are the same, 
and are obtained from the effective dielectric function e. The main contribution to the 
Casimir energy comes from frequencies in the vicinity of u ~ c/2L or wavelengths of the 
order of A ~ AirL. In order to calculate the Casimir force between composite materials the 
wavelength has to be larger than the typical size of the inhomogeneities a. That is, the 
separation between the slabs and the size of the inhomogeneities has to satisfy AttL > a. For 
example, if we consider a composite made of a homogeneous matrix and spherical inclusions 
of radius a = 20 nm, then the separations of the plates have to be larger than L > 1.5 nm. 

The Casimir force per unit area predicted by Eq. (9) calculated for the different EMA is 
shown in Figure (2). Using the common notation of the Casimir literature we have plotted 
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the reduction factor r), which is defined as the ratio of the force predicted using Eq. (9) 
to the force between perfect conductors (F = — ftc7r 2 /240L 4 ). The top panel of Figure 
(2) corresponds to a low filling fraction / = 0.015 where the Maxwell- Garnet, Bruggeman 
and lower Wiener bound give a similar result. The difference between Bruggeman model 
and the arithmetic average (upper Wiener bound) is of the order of 8% when the plates 
are 100 nm apart and 20% for plate separation of 300 nm. For a higher filling fraction of 
/ = 0.25 (bottom panel) the different models predict different values of the reduction factor, 
in particular at large separations. Also, for higher filling fractions a bigger magnitude of 
the reduction factor is observed. This is expected since the composite starts behaving more 
like a metal due to the higher fraction of Au. Also, we see that the difference between the 
Bruggeman and arithmetic average results increases, being of the order of 20% at 100 nm 
and up to 62% at 300 nm. These two models are singled out since they have been used in 
Casimir force calculations. 

The difference between the different models as a function of filling fraction can also be 
seen in Fig. (3), where the reduction factor is plotted as a function of filling fraction. The 
separation between the plates is kept fixed at L — lOOnm. For low filling fractions, the 
Bruggeman and Maxwell- Garnett models give the same results. The difference between 
them as the filling fraction increases is due to the insulator-metal transition predicted by 
the Bruggeman model. However, the difference with the arithmetic mean is significant for 
all filling fractions. 

IV. CONCLUSIONS 

The use of composite media to modify the Casimir force requires a careful study of the 
effective dielectric response of the system being used. As shown in this paper, the value 
of the force is very sensitive to the choice of the effective medium model. In the Casimir 
literature the upper Wiener bound and the Bruggeman approach are commonly used, but as 
we have shown, they can lead to very different theoretical predictions for the Casimir force. 
To make precise comparison between theory and experiments, in particular in high precision 
experiments, the choice of model needs to be carefully considered or misleading conclusions 
might arise. Before any measurement or calculation of the Casimir force is done, optical 
experiments are needed to determine which is the most suitable effective medium model for 
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a particular set up. 

V. APPENDIX. ELLIPSOIDAL INCLUSIONS 

If the inclusions are ellipsoids the Maxwell- Garnett and Bruggeman theories can be gen- 
eralized in the following way. The polarizability tensor of the ellipsoidal particle can be 
written as 



where a, b and c are the lengths of the ellipsoid semi-axes and Li are the depolarization 
factors. For a sphere L{ = 1/3 and for a prolate ellipsoid with a > b = c we have 



where e is the eccentricity of the elllipsoid. Replacing the previous expression in Eq.(3), the 
Maxwell-Garnett effective dielectric function for ellipsoids is obtained. 
Similarly, the Bruggeman formula can be generalized to ellipsoids as 
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FIG. 1: (Color online) Effective dielectric function plotted as a function of the filling fraction 
for a composite made of a S1O2 host with spherical Au inclusions. In the top panel the effective 
dielectric function is evaluated at the imaginary frequency ( = 0.02u p and at ( = 0.5lo p for the 
bottom graph. 
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FIG. 2: (Color online) Casimir force (normalized to the ideal case) as a function of plate separation 
between two composite plates made of a homogeneous matrix of S1O2 with spherical inclusions 
of radius a. The results in this plot are valid provided that condition AttL > a is satisfied. The 
top figure is for low filling fractions / = 0.015. The Maxwell-Garnett and Bruggeman models give 
results similar to the lower Wiener bound. For higher filling fraction (/ = 0.25) the difference in 
the force between the different models is much bigger in particular at large separations. 



14 



0.38 
0.36 
0.34 
0.32 
0.3 
0.28 
0.26 
0.24 



A dl III II 1 ICUO 

+ Bruggeman 






* 

* 


o Maxwell-Garnett 




* 

* 


+ 










* 




+ 


* 

* 




+ 

+ 


o 




* 




o 






+ 


o 




* + 

o 

+ 


o 






o 








* 6 




L= 


100 nm 



0.1 0.2 0.3 

filling fraction 



0.4 



0.5 



FIG. 3: (Color online) Reduction factor as a function of filling fraction when the two plates are at 
a fixed separation L = lOOnm. At low filling fraction the difference between the Bruggeman and 
Maxwell-Garnet models is negligible but as the amount of Au particles increases, the Bruggeman 
model correctly predict a more metallic behavior, increasing the force between the slabs. 
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